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Abstract 

In the presence of gratuitous inducers, the lac operon of Escherichia coli exhibits 
bistability. Most models in the literature assume that the inducer enters the cell 
via the carrier (permease), and exits by a diffusion-like process. The diffusive influx 
and carrier efflux are neglected. However, analysis of the data shows that in non- 
induced cells, the diffusive influx is comparable to the carrier influx, and in induced 
cells, the carrier efflux is 7 times the diffusive efflux. Since bistability entails the 
coexistence of steady states corresponding to both non- induced and induced cells, 
neither one of these fluxes can be ignored. Here, we formulate a model accounting for 
both fluxes. We show that: (a) The thresholds of bistability are profoundly affected 
by both fluxes. The diffusive influx reduces the on threshold by enhancing inducer 
accumulation in non-induced cells. The carrier efflux increases the off threshold by 
decreasing inducer accumulation in induced cells, (b) Simulations of the model with 
experimentally measured parameter values are in good agreement with the data for 
IPTG. However, there are discrepancies with respect to the data for TMG. They 
are most likely due to two features missing from the model, namely, the variation 
of the inducer exclusion effect and the specific growth rate with the lactose enzyme 
levels, (c) The steady states and thresholds obtained in the presence of both fluxes 
are well approximated by simple analytical expressions. These expressions provide a 
rigorous framework for the preliminary design of the lac genetic switch in synthetic 
biology. 
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1 Introduction 



The lac operon is a paradigm of the mechanisms controlling gene regulation. 
This interest was stimulated by the hope that insights into the mechanism 
of lac induction would shed light on the central problem of development, 
namely, the mechan i sm by which genetically identical cells acquire distinct 



pheno types (jMonodl . 119471 ; ISpiegelmanl . Il948l ) . 



Many of the early studies of the lac operon were concerned with the kinetics 
of enzyme induction in the presence of gratuitous inducers, such as thiomethyl 
galactoside (TMG) and isopropyl thiogalactoside (IPTG). The use of such 
inducers was important because they enabled the kinetics of enzyme synthesis 
to be separated from the masking effects of dilution. This was achieved by 
growing the cells in a medium containing a gratuitous inducer (which promotes 
enzyme synthesis, but not growth), and non-galactosidic carbon sources, such 
as succinate, glycerol, and glucose (which support growth, but not enzyme 
synthesis) . 

These early studies showed that the enzyme synthesis rate was not uniquely 
determined by the composition of the medium. If glucose and TMG were added 
simultaneously to a culture of E. coli ML30 growing on succinate, there was 
almost no synthesis of /3-galactosidase; however, if glucose was added 20 mins 
after the addition of TMG, /3-galactosidase was synthesized for up to 130 
generations (Fig. QJi). Thus, enzyme synthesis is bistable: Pre- induced cells 
remain induced, and non-induced cells remain non-induced. 

Cohn and Horibata found that the existence of bistability depended crucially 
upon the existence of lac pe rmease (LacY): It d i sappe ared in lacY~ (cryp- 



tic), but not lacL , mutants (ICohn and Horibatal . Il959l ). They proposed that 



bistability occurred due to the destabilizing effect of the positive feedback 
generated by lac permease. More precisely, they observed that: 



This permease not only accumulates galactosides intracellularly, but its syn- 
thesis is induced by the galactosides, which it accumulates. The permease 
may therefore be described as a "self-inducing" system, for which a general- 
ized schema might be drawn. 

This generalized schema, shown in Fig. [lb, was postulated before the molecular 
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Figure 1. (a) Bistability during g rowth of E. coli ML30 on succinate, glucose, and 
TMG (|Cohn and Horibatal . 11959 . Fig. 4). If glucose and TMG are added simultane- 
ously to a culture growing on succinate, there is no /3-galactosidase synthesis (♦). If 
glucose is added 20 min after the addition of TMG, the enzyme is synthesized (■) 
at a rate that is 50% of the rate observed when only T MG is added to the culture 
(a), (b) Kinetic scheme proposed by Cohn and Horibata ( Cohn and Horibata . 19591 . 
p. 611). Here, the enzyme represents lac permease; the substrate and product refer 
to extracellular and intracellular TMG, respectively. 



mech anism of induction was discovered by Jacob and Monod (j Jacob and Monod . 
196lh . 



Some years after the discovery of the induction mechanism, Babloyantz and 
Sanglier formulated a mathematical model, which was based on the Cohn- 
Horibata scheme, but described enzy me induction in terms of the no w well 
established Jacob- Monod mechanism (IBabloyantz and Sanglierl . Il972l ). They 
showed that the model yielded the bistability observed in experiments. Chung 
and Stephanopoulos formulated a similar model, the main difference being 
that repressor-ope rator and repressor-inducer bindin g were assumed to be in 
quasi-equilibrium (jChung and Stephanopoulosl . I1996T ). This model is given by 
the equations 



dx - t/ s - i ( 

^ ' s 1 x 1 ' s — ' s j£ | „ J x — x \ 



de 



r g e, r e 



1 + Klx 2 



K e 



(1) 
(2) 



l + a + Klx 2 ' c 

where x and s denote the intracellular and extracellular TMG concentrations, 
respectively; e denotes the permease activity; r g is the specific growth rate on 
the non-galactosidic carbon source(s); r s is the specific rate of TMG uptake by 
the permease; r x is the net specific rate of TMG efflux by diffusion; and r+, r~ 
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Figure 2. Bistability of the lac o peron during growth of E. coli K12 MG1655 on 
succinate and succinate + glucose (jOzbudak et all 12004 . Figs. 2b and 2c). (a) Bista- 
bility during growth on succinate and various concentrations of extracellular TMG. 
The (normalized) green fluorescence provides a measure of the steady state activ- 
ity of the lac operon. The upper (resp., lower) panel shows the green fluorescence 
observed when an induced (resp., non-induced) inoculum is grown exponentially on 
a mixture of succinate and various concentrations of extracellular TMG. Bistabil- 
ity occurs at extracellular TMG concentrations between the off and on thresholds 
of 3 and 30 fiM, respectively, defined as the extracellular TMG concentrations at 
which <5% of the cells are in their initial state, (b) Bistability persists even if glu- 
cose is added to the mixture of succinate and TMG. The lower (a) and upper (a) 
curves show the off and on thresholds observed at various glucose concentrations. 
The dashed curves show the thresholds predicted by the model with the parameter 
values in Table [TJ 

denote the specific rates of permease synthesis and degradation, respectively. 
The expression for r+ is based on a molecular model which assumes that the 
lac operon contains one o perator, and the lac r epressor contains two identi- 
cal inducer-binding sites (jYagil and Yagill . Il97ll ). The parameter, K x , is the 



association constant for the repressor-inducer binding; and a is jointly pro- 
portional to the intracellular repressor level and the association constant for 
repressor-operator binding. Evidently, a is a measure of the repression, defined 
as the ratio, r+l^^/ r+\ x=0 . 

Although the experiments done by Cohn and coworkers provided clear evi- 
dence of bistability during growth on TMG and succinate/glucose, they did 
not investigate the variation of the enzyme levels with TMG and glucose levels, 
and the effect of glucose on the key molecular mechanisms of lac regulation, 
namely, CRP-cAMP activation and inducer exclusion. Recently, Ozbudak et 
al. addressed these issues by inserting t wo reporter operons i nto the chromo- 
some of Escherichia coli K12 MG1655 flOzbudak et all 12004] ). 



The reporter lac operon, placed under the control of the lac promoter, coded 
for green fluorescent protein (GFP) instead of the lac enzymes. Thus, the 
GFP intensity of a cell provided a measure of the lac enzyme level. 
The reporter gat operon, placed under the control of the constitutive gat 
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Figure 3. Variation of the scaled red and green fluorescence intensities during ex- 
ponential growth of E. coli K12 MG1655 in the presence of succinate, TMG, and 
various concentrations of glucose, (a) The scaled r ed and green fluoresce nce inten- 



sities are identical at all glucose concentrations (jOzbudak et all 12004 Fig. 3b). 



(b) The scaled red fluoresc ence intensity decreases with the concentration of glu- 



cose (jOzbudak et all [200J, Fig. 3a) . The curve shows the fit to Eq. ([25 



promoter, coded for the red fluorescent protein (RFP) instead of the gat 
enzymes. The RFP intensity of a cell provided a measure of its CRP-cAMP 
level. 



To assess the effect of various TMG and glucose levels on lac bistability, they 
exposed non-induced and induced cells to a fixed concentration of succinate, 
and various concentrations of TMG and glucose. It was observed that: 

(1) When the cells were grown in the presence of succinate and various con- 
centrations of TMG, they exhibited bistability (Fig. [2k). 

(2) This bistability persisted even if glucose was added to the mixture of 
succinate and TMG, but the thresholds increased with the concentration 
of extracellular glucose (Fig. Eb). 

They also showed that both observations were mirrored by the bifurcation di- 
agram for the modified Chung-Stephanopoulos model, obtained by neglecting 
the diffusive influx (i.e., the term, k x s, in Eq. (QQ) was assumed to be zero). 

To quantify the lac repression due to glucose-mediated reduction of CRP- 
cAMP levels, they measured the average green and red fluorescence intensi- 
ties of cells grown overnight in a medium containing a high concentration of 
TMG, a fixed concentration of succinate, and various concentrations of glu- 
cose. They found that the green and red fluorescence intensities, scaled by the 
corresponding values observed in a medium containing no glucose, were equal 
at all glucose concentrations (Fig. [3k). Furthermore, the normalized RFP in- 
tensity decreased 5-fold at saturating glucose concentrations (Fig. \3jp). Ozbu- 
dak et al. assumed that the GFP intensity is proportional to the maximum 
lac induction rate, and the RFP intensity is proportional to the CRP-cAMP 
level. Given these two assumptions, it was concluded from Figs. [3k and b 
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Figure 4. The kinetics of induction and diffusive influx in the cryptic mutant, E. coli 
ML3 (jHerzenbergl . Il959l . Figs. 1 and 4). (a) Variation of the scaled induction rate 
with extracellular TMG and IPTG levels. The full curves show the fits to the data 
obtained with Eq. ([2D and the parameter values a = 40, a = 1200, K~ l = 24 fj,M 
for IPTG, and K' 1 = 470 /xM for TMG. (b) Existence and kinetics of the diffusive 
influx. The cryptic mutant hydrolyzes ONPG, and the hydrolysis rate increases 
linearly with the extracellular TMG level. The permease activity is expressed in 
fimols min^ 1 gdw -1 . 

that the maximum lac induction rate is proportional the CRP-cAMP level, 
and decreases 5-fold at saturating concentrations of glucose. This is somewhat 
higher than the maximum possible cAMP-mediated repression in E. coli K12 
MG1655. Recently, Kuhlman et al. found that during growth of this strain on 
saturating concentrations of 0.5% glucose and 1 mM IPTG, variation of the 
extracellular cAMP level from to 10 mM pr oduces no more than a 3-fold 
change in lac expression (IKuhlman et all 120071 . Fig. IB). 



The modified Chung-Stephanopoulos model considered by Ozbudak et al. pro- 
vides an excellent point of departure for quantifying lac bistability. However, 
it suffers from three deficiencies. 



The induction kinetics are incorrect The Yagil and Yagil model of lac 
induction is not consistent with the structure of the lac operon and repres- 
sor. Indeed, the lac operon contains two auxiliary operators, O2 and O3, in 
addition to the main opera tor, O i , and the lac repressor contains four (not 
two) inducer-binding sites (jLewisl . 120051 ). Furthermore, these structural fea- 
tures play a cruc ial role in the formation of DNA loops, t he key determinant s 
of lac repression dOehler et allll990lll994h and induction (lOehler et all 120061 ). 
Molecular models taking due account of the 3 operators and 4 inducer-binding 
sites yield the lac induction rate 



K 



1 + a/ (1 + K x xf + d/(l + K x xf ' 



(3) 



where K x is the association constant for repressor-inducer binding, and a, a 
are related to the the repression stemming from repressor-operator binding and 
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DNA looping, respectively (IKuhlman et all 120071 ; INarangl . 120071 ). Equ ation ([31 ) 
yield s good fits to the induction curves obtained with various strains (INarangl . 
20071 . Fig. 11) and gratuitous inducers (Fig. Hk). 



The diffusive influx is absent It is well known that there is a diffusive 
influx. Herzenberg showed the existence of the diffusive influx by measuring 
the ONPG hydrolysis rate in wild-type and cryptic cellsJZI He found that even 
cryptic cells hydrolyzed ONPG, and the hydrolysis rate increased linearly with 
the extracellular ONPG concentration (Fig. Hb)- It follows that gratuitous 
inducers can enter the cells by a permease-independent mechanism, which can 
formally be described by diffusive (first-order) kinetics. 

Now, it may be permissible to neglect the diffusive influx in induced cells 
because such cells contain such large permease levels that the inducer is im- 
ported almost entirely by the permease. Ho wever, non-induced cells con tain 



only ~0.1% of the permease in induced cells (IMaloney and Wilson! . Il973l ). In 
such cells, it seems unlikely that the diffusive influx is negligible. As we show 
below, this conjecture is supported by the data. 



The carrier efflux is absent The existence of the carrie r efflux was sug- 
gested in early studies with gratuitous inducers (IKochl . 19641) . However, Mal- 



oney and Wilson were the first to quantify its effect (IMaloney and Wilson! 
19731 ). To this end, they measured the steady state intracellular TMG levels 
in cells induced to various levels, and then exposed to chloramphenicol plus 
0.24 niM extracellular TMG. The data obtained, shown in Fig. [5^, cannot 
be reconciled with models, such as the Chung-Stephanopoulos model, that 
account only for the carrier influx and diffusion. Indeed, since enzyme synthe- 
sis is blocked in the presence of chloramphenicol, the Chung-Stephanopoulos 
model implies that the steady state intracellular TMG level is given by the 
expression 



x 



I — 

- krr. 



V = V*e, 



(4) 



K t + s 

obtained from ([I]) by ignoring the negligibly small dilution term. It follows 
that if cells containing various permease levels are exposed to chloramphenicol 
and a fixed extracellular TMG level, the intracellular TMG level, x, should 
increase linearly with V, the specific activity of the permease. However, the 
experiments show that x increases linearly only if V is small — it saturates 
at large values of V (closed circles in Fig. [5^). 



2 The ONPG hydrolysis rate in int act cells is proportional to the rate at which 
ONPG (jHerzenbergl . Il959l ) or TMG (jMalonev and Wilsonl . Il973l . Fig. 1) permeate 
the cell. 
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Figure 5. Existence and quantification of the carrier efflux, (a) When E. coli K12 
CA8000, induced to various levels, is exposed to chloramphenicol and 0.24 M extra- 
cellular TMG, the ste ady state intracellular TMG level increases nonlinear ly with 
the permease activity ( Malonev and Wilson . 19731 . Fig. 5). The permease activity is 



normalized by the activity of fully induced cells (43.5 //moles TMG per min per mL 
of cell water) . The full curve shows the fit obtained from Eq. (0) with the experimen- 
tally measured parameters, s = 0.24 mM, K\ = 0.8 mM, k x = 0.14 min -1 , and the 
best fit parameter estimate, K 2 = 84 mM. The dashed line shows the intracellular 
TMG level predicted by Eq. (T4]). (b) Flux ratios estimated from ([5]). In fully induced 
cells, the carrier efflux is ~ 7 times the diffusive efflux. In non-induced cells, which 
have a scaled permease activity of ~0.001, the diffusive influx is comparable to the 
carrier influx. 

To resolve this discrepancy, Maloney and Wilson postulated that the permease 
catalyzes both influx and efflux, i.e., the steady state intracellular TMG level 
satisfies the equation 

° = V (l^-l^—)- k *( x - s )> ( 5 ) 
\Ki + s K 2 + xJ 

where the term, Vx/(K 2 + x), accounts for the carrier efflux. They measured 
K±, k x , and showed that ([5]) yields a good fit to the data, provided K 2 is chosen 
to be 84 mM (full curve in Fig. [5^). 

With the help of Eq. ([5]), they explained their data as follows. If the cells 
have low enzyme levels, the intracellular TMG levels are so small compared 
to K 2 that the carrier efflux is negligible. Thus, ([5]) reduces to (@J, and the 
intracellular level increases linearly with V (dashed line in Fig. [5^). However, 
at large enzyme levels, the intracellular TMG level becomes so high that the 
steady state is completely determined by the first two terms of ([5]), namely, the 
carrier influx and efflux. The intracellular TMG level is therefore independent 
of the enzyme activity. 

The foregoing physical argument assumes that the carrier efflux is negligible 
in non-induced cells, and diffusion is negligible in induced cells. These assump- 
tions are consistent with the mathematical model. Indeed, given Eq. ([5]) and 
the above parameter values, we can plot the two ratios, carrier efflux:diffusive 
efflux and diffusive influxxarrier influx, as a function of the enzyme activity 
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(Fig. [5b). The plots show that in induced cells, the carrier influx is ~1000 
times the diffusive influx (full line), and the carrier efflux is almost ~7 times 
the diffusive efflux (dashed line). In non- induced cells, on the other hand, the 
carrier efflux is negligible compared to the diffusive efflux, but the diffusive 
influx is comparable to the carrier influx. 

It is therefore clear that the carrier efflux cannot be ignored in induced cells, 
and the diffusive influx cannot be neglected in non-induced cells. Since bista- 
bility entails the coexistence of steady states corresponding to both induced 
and non-induced cells, neither the carrier efflux nor the diffusive influx can 
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The goal of this work is analyze this effect. 

Recently, we studied the effect of DN A looping on the dynamics o f the mod- 
ified Chung-Stephanoupoulos model (INarang and Pilyuginl . 120071 ). Here, we 
formulate and analyze an extended model of lac transcription that takes due 
account of DNA looping as well as the diffusive influx and the carrier efflux. 
We find that: 



(1) The diffusive influx has no effect on the off threshold, but it significantly 
reduces the on threshold. 

(2) The carrier efflux has no effect on the on threshold, but it significantly 
increases the off threshold. 



We show that the on and off thresholds can be represented by simple analytical 
expressions. 



2 Model 



Fig. [6] shows the kinetic scheme of the model. Here, E denotes lac permease; 
S and X denote the extracellular and intracellular inducer, respectively; GFP 
and RFP denote the green and red fluorescent protein synthesized by the 
reporter lac and constitutive gat operons, respectively; and GLU, SUC, GLY 
denote glucose, succinate, and glycerol. We assume that: 

(1) The concentrations of the permease, intracellular inducer, GFP, and RFP, 
denoted e, x, g, and r, respectively, are based on the volume of the cells 
(mol L _1 ). The concentrations of the extracellular inducer and glucose, 
denoted s and G, respectively, are based on the volume of the reactor 
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Figure 6. Kinetic scheme of the model. 

(mol L^ 1 ). The rates of all the processes are based on the volume of the 
cells (mol h _1 L _1 ). We shall use the term specific to emphasize this 
point. 

The choice of these units implies that if the concentration of any intra- 
cellular component, Z, is z mol L _1 , then the evolution of z in batch 
cultures is given by 

dz , 

dt z z 9 

where r+ and r~ denote the specific rates of synthesis and degradation 
of Z in mol h -1 L -1 , and r g is the specific growth rate in h 
(2) The specific growth rate is given by 



r g = r g0 (f)g(G), 

where r g o denotes the specific growth rate in the absence of glucose (i.e., 
in the presence of pure succinate), and 4> g {G) is an increasing function 
of the extracellular glucose concentration, G, satisfying (f) g (0) = 1, which 
accounts for the increase of the specific growth rate produced by addition 
of glucose to a culture growing on succinate. 

The saturation constant for carrier efflux, is so large that xj (K2+X) m 
x/K 2 . Hence, the net carrier-mediated specific uptake rate of S is 



r s = V s0 4> s (G)e 



x 



where V s oe is the specific activity of the permease in the absence of glu- 
cose, and 4> S (G) is a decreasing function of G such that <f> s (0) = 1, which 
accounts for the reduction of the specific activity due to inducer exclusion. 
The net specific rate of expulsion of X by diffusion follows the kinetics 



r x — k x (x s) . 
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(5) The intracellular inducer stimulates the transcription of the native and 
reporter lac operons, resulting in the synthesis of the lac enzymes and 
GFP, respectively. 

(a) The specific synthesis rate of the permease, E, is 



T e = V e Q<j> e ,lac{G) 



1 + a/ (1 + K x x) +&/(! + K x x 



A ; 



(6) 



where V e o is the specific synthesis rate in fully induced cells growing 
on succinate (G — 0); 4> e ,iac{G) is a decreasing function of G satisfying 
e ,Zac(O) = 1, which accounts for the repression of transcription due 
to reduction of CRP-cAMP levels in the presence of glucose; K x is 
the association constant for repressor- inducer binding; and a, a char- 
acterize the repression due to repressor-operator binding and DNA 
looping, respectively. 

The values of a and a are 20-50 and ~1200, respectively, which reflect 
the fact that more than 95% of t he total repression, 1 + a + a, is due 
to the formation of DNA loops flOehler et all . ll990L Il994j ). The dis- 
sociati on constant for repressor-inducer binding, K~ l , is 7-30 /iM for 
IPTG flNarand . I2007L Table 1). It is 10 times higher f or TMG, based 



on a simple analysis of the data with Scatchard plots (IBarkley et al. 



19751 . Table I). A similar estimate is obtained when Eq. ([3]) is used 
to fit the data shown in Fig. H^. Assuming a = 40 and a = 1200, the 
best fits to this data are obtained when the dissociation constants for 
IPTG and TMG are assumed to be 24 and 470 /xM, respectively, 
(b) The specific synthesis rate of GFP is 



r GFP = V e0 (f) e: i ac (G) 



1 



1 + a/ (1 + K x x) + & G / (1 + K x x 



A - 



where «g < ot. This expression is obtained by assuming that the 
promoters of the reporter and native lac operons are identical, so 
that both operons have the same maximum specific synthesis rates 
in the absence of glucose, and are subject to the very regulation 
by CRP-cAMP and repressor-operator binding. However, the two 
operons differ with respect to regulation by DNA looping because the 
reporter lac operon lacks the auxiliary operator, 2 . This precludes 
the formation of DNA loops by interaction between 0\ and O2. 
Ozbudak et al. estimate d the total rep r ession of the lac reporter, 
1 + a + c\g, to be 170 (jOzbudak et all 120041 ) . which implies that 
a G ss 130. 

(6) The specific synthesis rate of RFP follows the constitutive kinetics 

r RFP = VrO(fie,gat(G), 

where V r0 is the specific synthesis rate of RFP in the absence of glucose, 
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and 4> e ,gat(G) is a decreasing function of G such that e , 9 at(O) = 1, which 
accounts for the CRP-cAMP effect exerted in the presence of glucose. 
(7) Protein degradation is negligible. 

Given these assumptions, the mass balances yield 

= V s0 (f) s e ( — — -k x (x-s)- r g0 (j)gX, (7) 



dt \Kt + s K 2 , 

de 1 

* = V ^-" 1 + a/(1 + K lX f + a/( 1 + K lX f - (8) 

it = v ^-^\ +a/{1+ K J + &ol (1 + KiX f - r ^ < 9 > 

dv 

-j~ t = Ko0e,got - r g^g r - (10) 

It is convenient to rescale these equations. Ozbudak et al. scaled the GFP 
and RFP intensities by the corresponding intensities observed during steady 
exponential growth on succinate and excess TMG, but no glucose. Under these 
conditions, Eqs. (!8l)-(fT0l) imply that e, g, and r have the values V e o/r g0 , V e o/r g0 , 
and V r o/r g0 , respectively. Thus, we are led to define the dimensionless variables 

e g _ r 

Veo/rgo 1 7 V e0 /r g o ,P V r0 /r g0 ' 

If we scale the remaining variables, x and t, as follows, 

x t 



X 



fC-l> — r -i' 



we arrive at the dimensionless equations 

r x ^- = 5 m0 (f) s e ( — — ) - (x ~ cr ) - T x <p g x (11) 

dr \Kx + a K2/ 

d/€ X 

-r = <Pe,iaJ{x) ~ 9 e, fix) = Z JT, 72 rrr, 74' ( 12 ) 

dr 1 + «/ (1 + x) + «/ (1 + x) 

d<v 1 

dr 1 + «/ (1 H- x) +Og/(1+x) 



dp _ , 
dr 9e ' 9at ' 

with the dimensionless parameters, 



(14) 



r g o s Ki V s0 (V e0 /r g0 ) /k x 
T ^T< a =K?-*'= K?' *- s • (15) 



X X X 



Here, r x is the ratio of the time constant for inducer expulsion by diffusion 
(/c" 1 ) relative to the time constant for growth on succinate (r^, 1 ); a is the ex- 
tracellular inducer concentration, measured in units of K~ x \ K\ and k 2 are the 
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saturation constants for carrier uptake and expulsion, respectively, measured 
in units of K~ x 



X 



The parameter, 5 m0 , is related to the steady state intracellular inducer level 
attained in fully induced cells exposed to chloramphenicol, succinate, and sat- 
urating inducer levels, if the carrier efflux is somehow abolished. Indeed, since 
r 9j0 = 0.4 h _1 , k x = 0.14 min -1 , and (fi g < 2, the dilution term, T x (p g x, in (JTT1) 
is negligibly small compared to the diffusive efflux, \. The steady states are 
therefore given by the relations 

_ cr + <5 m o0 s eq-/(/ti + a) 
1 + o m0 <p s e/K 2 

e=% 2£ /(x), (17) 
<Pg 

7 = %^MX), (18) 

p =^fL. (19) 



Equation fTT6l) provides the steady state intracellular inducer level in cells 
exposed to fixed concentrations of chloramphenicol and succinate, and any 
given concentrations of extracellular inducer and glucose. If the carrier efflux 
is somehow abolished (k 2 — > oo), and the medium contains no glucose (<p s = 1), 
the steady state intracellular inducer level is given by the equation 

X = v + 8 m0 e—?— . (20) 
k 1 + a 

If the cells are fully induced (e = 1) and the extracellular inducer level is 
saturating (cr ^> ki), the above equation becomes 



X = cr + 5m0 ~ <W 

where the first term was neglected because in induced cells, the diffusive influx 
is vanishingly small compared to the carrier influx. Thus, 5 m o is the steady 
state intracellular inducer level, measured in units of K' 1 , that is attained in 
fully induced cells exposed to chloramphenicol, succinate, and saturating levels 
of extracellular inducer, provided the carrier efflux is somehow abolished. 

We can estimate the value of 8 m oK~ l by observing that (1201) is the dimension- 
less analog of (jl]). Indeed, multiplying (120]) by K~ l yields the relation 

x = s + K^S^e 
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which, upon comparison with (j3J), gives 



fc-lA - V 
J^x f m o e — ~T~ ■ 

Since V = 43.5 /imoles min -1 mL _1 in fully induced cells, we conclude that 
dmoK- 1 = V/k x = 310 mM. 

If the medium contains the inducer, succinate, and possibly glucose, but no 
chloramphenicol, enzyme synthesis is not blocked. Instead, it responds to the 
prevailing intracellular inducer concentration, and the enzyme level evolves to- 
ward the steady state given by (ITT)) . In this case, the steady state intracellular 
inducer level satisfies the relation 

cr + 5 m (G)f(x)(r/(K 1 + (x) _ _ (p s ^i ac 

X = 7— ; (n , f( w , o m (G) = <W(G), <j){G) = — — , 21 

l + d m {G)f{x)/K 2 <p g 

obtained by substituting ( ITT)) in ( fl6l) . Here, 4>{G) < 1 captures the cumulative 
effect of glucose due to inducer exclusion, catabolite repression, and enzyme 
dilution; and S m (G) is a measure of the steady state intracellular inducer level 
in fully induced cells exposed to succinate, glucose, and saturating inducer 
levels. 

We are particularly interested in: 

(1) The variation of the steady states with the extracellular inducer concen- 
tration at any given glucose level (for example, G = 0, in Fig. [2^). This 
is completely determined by Eq. (121)) : Given any a and G, the equation 
can be solved for %, which can then be substituted in (|TT)) - (IT5)) to ob- 
tain the corresponding e and 7. We shall refer to ()2ip as the equilibrium 
condition. 

(2) The variation of the thresholds with the extracellular glucose concen- 
tration (Fig. |2b). These thresholds are the points at which the steady 
state bifurcates (i.e., changes its stability or multiplicity). It is shown in 
Appendix [X] that a steady state bifurcates only if it satisfies the relation 



— if + Xfx) - ~^—fx 

K>2 K\ T <J 



(22) 



We shall refer to this equation as the bifurcation condition. 



Thus, the steady states are completely determined by the equilibrium con- 
dition, whereas the thresholds must satisfy the equilibrium and bifurcation 
conditions. 
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3 Results and Discussion 



3. 1 The effect of dilution 



To simulate the effect of the diffusion influx and the carrier efflux, we need the 
function, 4>{G) = ^ a </> e /<j> g . In the course of estimating this function, we shall 
also resolve the discrepancy between the magnitudes of the cAMP-mediated 
repression reported by Ozbudak et al. and Kuhlman et al. 

It follows from f[T8l) and ffT9|) that when the cells are grown in the presence of 
excess TMG, the ratio of the steady state GFP and RFP intensities is given 
by 

7 _ 4>e,lac(G) 
P 4>e,gat(G) 

Since this ratio was observed to be 1 all glucose concentrations (Fig. [3k), the 
lac and gat operons respond identically to CRP-cAMP, i.e., 

4>e,lac = 4>e,gat = <j>e, Say. (23) 

Moreover, the steady state RFP and GFP intensities are given by the equation 

P = l = ^. (24) 

It follows that the 5-fold decline of p and 7 at saturating glucose concentra- 
tions, shown in Fig. [3b, represents the combined effect of reduced CRP-cAMP 
levels and enhanced dilution (as opposed to the the sole effect of reduced 
CRP-cAMP levels). 

The effect of the reduced CRP-cAMP level, <j) e , cannot be estimated unless 
the effect of enhanced dilution, <f) g , is known. Ozbudak et al. did not report 
the specific growth rates at va rious glucose concentrations, but experiments 



show that (INarang et all Il997l ) : 



• When E. coli K12 is grown on saturating levels of succinate + glucose, the 
cells consume only glucose during the first exponential growth phase. 

• The maximum specific growth rate on glucose (0.74 h -1 ) is roughly 2 times 
that on succinate (0.44 h^ 1 ). 

It follows that at saturating glucose concentrations, <p g ~ 2. Since p ~ 1/5 
under these conditions, (124"}) yields e = pcp g = 2/5. Thus, roughly half of the 5- 
fold decline in Fig. [3b is due to the enhanced dilution rate at saturating glucose 
concentrations. Reduction of the CRP-cAMP levels accounts for the remaining 
2.5-fold decline, which is consistent with the data obtained by Kuhlman et al. 
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Table 1 

Parameter values used in the simulations. All the parameter values, except K~ , 
are assumed to be the same for TMG and IPTG. 



Parameter 


Value and reference 


Kl 


0.68 mM COzbudak et al.. 2004) 


K 2 


84 mM (Malonev and Wilson. 1973) 


K? 


0.3 mM for TMG fHerzenbertr. 1959) 




0.007 mM for IPTG COehler et al.. 2006) 


a 


40 fOehler et al.. 1994) 


a 


1200 fOehler et al.. 1994) 


a G 


130 rOzbudak et al.. 2004) 


4*e/<f>g 


See Eq. ([25]) 


4>s 


See Eq. ((26D 




310 IKZ 1 fMalonev and Wilson. 1973) 



In view of (1241) . we can assume that 

^ = 1 - 0.84- 



G 



7.6 + G" 

which represents the best fit obtained to the data in Fig. [3b. 



(25) 



It remains to specify the function, <f) s , which characterizes the intensity of in- 
ducer exclusion. We assume that the saturation constants for inducer exclusion 
and cAMP activation/dilution are the same (7.6 /iM). Experiments show that 
the lacto se uptake rate decrease s 2-fol d in the presence of high glucose concen- 
trations (jMcGinnis and Paigenl . 1 19691 . Fig. 3). Thus, we are led to postulate 
the expression, 

We shall discuss the implications of the foregoing assumption later on. 



3.2 The effect of the diffusive influx and carrier efflux 



All the parameters required to study the effect of the diffusive influx and carrier 
efflux are now available (Tabled]). Most of the simulations shown below were 
done with the parameter values for TMG. However, at the end of this section, 
we shall show the simulation for IPTG, and compare it with the experimental 
data. 
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(a) (b) 

Figure 7. Variation of the steady states with the extracellular TMG level, a, at 
G = in the absence of diffusive influx and carrier efflux, (a) Enzyme activity (full 
line) and green fluorescence intensity (dashed line), (b) Intracellular TMG level. 

We shall begin by considering the base case of no diffusive influx or carrier 
efflux. To this base case, we shall add the diffusive influx and the carrier efflux, 
one at a time, to examine their effects on the steady states and the thresholds. 
Finally, we shall show that that the general model containing both fluxes is 
essentially a composite of the two special cases accounting for the diffusive 
influx and carrier efflux. 



3.2.1 No diffusive influx or carrier influx 

In the absence of the diffusive influx and the carrier efflux, the equilibrium 
condition (I21j) becomes 

X = S m (G)f( X )^—. (27) 
Ki + a 

We can determine the variation of the steady states with a at any given G by 
solving this equation for x, and substituting it in (fTTj) and ffl~8l . However, this 
approach can only be implemented numerically. It is more convenient to solve 
f l27|) for cr, thus obtaining 

a( X ,G) = K 1 - f . (28) 

o m {G)f{x) - x 

The variation of the steady state e, 7, and x with a is then given by the 
parametric curves, {cr(x, G), e(x, G)}, {a(x, G), t(x, G)}, and {a(x, G), x)}, 
respectively, where e(x, G), j(x,G), and cr(x, G) are given by f|T7|) . ffTSl . and 
(125]). respectively. 

Figure [7^ shows the variation of the steady state e and 7 with a at G = 0. 
Both variables are bistable over the very same ~10-fold range of extracellular 
TMG concentrations. Furthermore, the GFP intensity and enzyme activity of 
the induced cells are identical (e, 7 « 1), but the GFP intensity of non-induced 
cells is significantly higher than their enzyme activity. 
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(c) (d) 

Figure 8. Variation of the on and off thresholds in the absence of diffusive influx and 
carrier efflux, (a) Variation of the on threshold (upper full curve) and off threshold 
(lower full curve) with S m . The dashed vertical lines represent the values of 5 m at 
which the thresholds become infinitely large. The dotted vertical lines represent the 
values of S m at G = and G = 1000 fiM. (b) If 5 m = 10 < l// x (x 2 ), there is 
no bistability. (c) If 5 m = 80, which lies between l// x (xi) and l/f x (x2), there is 
bistability, but no on threshold, (d) Variation of the off threshold (lower curve) and 
on threshold (upper curve) with the glucose concentration, G. 

The foregoing trends can be understood in terms of the variation of the intra- 
cellular TMG levels with o (Fig. [Tb). In non- induced cells, the GFP intensity 
is significantly higher than the enzyme activity because the intracellular TMG 
levels are so small compared to K~ x that the reporter and native lac operons 
are transcribed at their basal rates. Since the reporter lac operon lacks O2, 
its basal transcription rate, 1/(1 + a + &c) = 1/170, is ~7 times higher than 
that of the native lac operon, 1/(1 + a + a) = 1/1241. In induced cells, on 
the other hand, the inducer levels are so high that the repressor is completely 
inactivated, and both operons are transcribed at the very same (maximal) 
rates. Finally, the thresholds for e and 7 are the same because f|T71) and f|T8l) 
imply that e and 7 are increasing functions of \. Hence 

do do d\ do do d\ 
de dx de ' d^ d\ d'y 

are zero precisely when do/d\ = 0. Since do/dx is zero at o w 0.03 and 
o « 0.3, so are do /de and do jd^. 
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We are particularly interested in the variation of the on and off thresholds with 
G, the concentration of extracellular glucose. In this regard, it is worth noting 
two points. First, since e and 7 always have the same thresholds, it suffices to 
focus on either one of these variables. In the simulations that follow, we present 
only the enzyme activities. Second, although we are interested in the variation 
of the thresholds with G, we shall begin by determining their variation with 
5 m . As we show below, this immediately yields their variation with G because 
5 m is a monotonically decreasing function of G. 

In the absence of the carrier efflux, the bifurcation condition (1221) becomes 



The thresholds satisfy (1271) and ([29]) . which imply that 

X ~ 7- = 0. (30) 
jx 

It is shown in Appendix [B] that (|30l) has exactly two positive roots, xi < X2, 
and they satisfy the relation, f x (xi) < fxiXz)- The loci of the on and off 
thresholds are therefore given by the equations 

<W x (xi)-i' <Wx<A2)-l' 
which define two decreasing curves on the <5 m cr-plane (Fig. [8^). Furthermore: 

• The on and off thresholds become infinitely large when 8 m approaches 
l// x (Xi) and l/f x (x2), respectively (dashed lines in Fig. (Hh). 

• When 5 m > l/f x {Xi) > V/xM, d3B> reduces to 



a 




/v(Xa)/ 



Thus, the curves are hyperbolic, and the ratio of the on threshold to the off 
threshold equals the constant, fxix?)/ fx(Xi)- 

Taken together, these results imply that if 5 m > 1/ f x (xi), there is bistability 
with finite on and off thresholds (Fig. |7J); if 5 m < l// x (x 2 ), there is no bista- 
bility (Fig. [8b); and finally, if l// x (% 2 ) < 5 m < l// x (xi), there is bistability 
with a finite off threshold, but no on threshold (Fig. [8b). We shall elaborate 
on the last case shortly. 

The variation of the off and on thresholds with the glucose level, G, follows 
immediately from Fig. [HK Indeed, as G increases from to 1000 /iM, 5 m = 
S m o4> decreases from 5 m0 to 5 moo = S m o(f>\ G=1000 , shown by the dotted vertical 
lines in Fig. [Hk. It follows that as G increases, the on and off thresholds as well 
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as the ratio of the on to off threshold increase. Furthermore, the off threshold 
exists for all < G < 1000 //M, whereas the on threshold ceases to exist for 
some G < 1000 fM. 

We can get explicit expressions for the variation of the thresholds with G by 
letting S m = 5 m0 (p(G) in fl3Tl) . thus obtaining 

O = t , , ^ „ ; ; -, O 



W(G)/ X (xi) - 1' W(G)/ X (x 2 ) - 1' 

Since decreases with G, a increases with G. Thus, on the Gcx-plane, the 
thresholds are represented by two increasing curves (Fig. [HJd) . As expected, 
these curves diverge with increasing glucose levels, with the on threshold be- 
coming unbounded at G ~ 100 fiM. 

The disappearance of the on t hreshold at intermediate values of 5 m is an exam- 



ple of irreversible bistability (ILaurent and Kellershohnl . Il999l . p. 421). Indeed 



it is evident from Fig. [8b that synthesis of the enzymes in induced cells can be 
switched off by gradually decreasing the extracellular inducer concentration. 
However, once the enzyme synthesis has been switched off in this manner, it 
cannot switched on again by any gradual increase of the extracellular inducer 
concentration. Irreversible bistability is thought to play an important role in 
development because it furnishes a mechanism for changing the cellular pheno- 
type permanently. However, it is biologically implausible in the case of the lac 
operon. Due to the existence of the diffusive influx, non- induced (and in fact, 
even cryptic cells) can always be induced by exposing them to sufficiently 
high extracellular inducer levels (Fig. H^). We show below that irreversible 
bistability disappears in the presence of the diffusive influx. 



3.2.2 Diffusive influx, but no carrier efflux 

In the absence of the carrier efflux, (|2"T|) becomes 

x = a + 5 m (G)f- 



Kx + O' 



which can be solved for a to obtain 



Ki + S m f - X) + V («1 + $mf ~ X) + ^KxX 



(32) 



(33) 



The variation of the steady state e, 7, and x with a is given by the parametric 
curves, {cr(x, G), e(%, G)} and {c(x, G),x)}, respectively, where e(x, G) and 
a(x,G) are given by ffPTl) and fl33|) . respectively. 

Figure [9^ compares the variation of the steady state permease activities with 
cr at G = in the presence (full curve) and absence (dashed curve) of the 
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(b) 



Figure 9. Comparison of the steady states at G = in the presence (full lines) 
and absence (dashed lines) of the diffusive influx, (a) Variation of the steady state 
permease activity, e, with the extracellular TMG level, a. (b) Variation of the steady 
state intracellular TMG concentration, x, with the extracellular TMG level, a. 

diffusive flux. It is evident that the diffusive flux has no effect on the off 
threshold, but it significantly reduces the on threshold. The physical reason 
for this is as follows. The off threshold is a property of the induced cells. 
These cells contain such high permease levels that the diffusive influx makes 
virtually no contribution to the accumulation of TMG (see upper branches 
of the curves in Fig. [9b). In contrast, the on threshold is a property of the 
non-induced cells, which contain such low enzyme levels that the diffusive 
influx significantly improves the accumulation of TMG (see lower branches of 
the curves in Fig. [9b). This enhanced accumulation of TMG decreases the on 
threshold. 

We have shown above that in the absence of glucose, the diffusive influx influ- 
ences only the on threshold. Is this true even in the presence of glucose? We 
can address this question by constructing the loci of the on and off thresholds 
on the # m <7-plane. To this end, observe that the thresholds satisfy the equilib- 
rium condition ([32]) and the bifurcation condition fl29l) . These two equations 
can be solved for a and 5 m to obtain the relations 



which provide the parametric representation, {d m (x)i a (x)}i °f the thresholds 
on the 5 m , a-plane. The full curve in Fig. [TUk shows a plot of this represen- 
tation. It consists of an upper branch and a lower branch, which represent 
the on and off thresholds, respectively. Since these two branches meet in a 
cusp, bistability is reversible whenever it exists: Either both or none of the 
thresholds occur at any given S m . As expected, irreversible bistability, which 
corresponds to the existence of only one of the thresholds, disappears in the 
presence of the diffusive flux. Importantly, this conclusion does not depend 



°{x) 




°{x) 



(34) 



(35) 



21 




Figure 10. Variation of the on and off thresholds in the presence of the diffusive 
influx, (a) Variation of the on and off thresholds with 5 m . The lower and upper full 
curves, defined by (|3"3j) and ([35]) . represent the off and on thresholds, respectively, 
in the presence of the diffusive influx. The full circle shows the cusp at which the 
thresholds merge. The lower and upper dashed curves, defined by (|3ip . represent 
the off and on thresholds, respectively, in the absence of the diffusive influx. The 
dotted vertical lines show the values of 5 m at G = and G = 1000 fiM. (b) Variation 
of the off threshold (lower curve) and on threshold (upper curve) with the glucose 
concentration, G. 



on the precise values of the parameter values chosen for the simulations. It 
is shown in Appendix O that the geometry of the curves representing the 
thresholds is preserved, regardless of the parameter values. 

To address the question posed above, it is useful to compare the thresholds in 
the presence of the diffusive influx with those obtained in its absence (dashed 
curves in Fig. [TUh). This comparison shows that there are 3 distinct regimes. 
When 5 m > 2000, the full and dashed curves coincide, i.e., the diffusive in- 
flux has no effect on both thresholds. At large 5 m , the permease is so active 
that even in the non-induced cells, the diffusive influx makes no contribution 
to inducer accumulation. When 100 < <5 m <2000, the lower full and dashed 
curves coincide, but the upper full curve is significantly lower than the upper 
dashed curve. The influence of the diffusive influx is therefore identical to that 
observed at G = 0: It has no effect on the off threshold, and significantly de- 
creases the on threshold. At these intermediate 5 m , the enzyme activity of the 
induced cells remains so large that the diffusive influx makes no contribution 
to inducer accumulation, but the enzyme activity of the non-induced cells be- 
comes so small that the diffusive influx does enhance inducer accumulation. 
Finally, when 5 m < 100, the enzyme activity becomes so small that the diffu- 
sive flux contributes significantly to inducer accumulation in both non-induced 
and induced cells. In fact, at sufficiently small 5 m , the inducer accumulates al- 
most entirely by diffusion, i.e., x ~ a - Since the intracellular TMG level is 
independent of the enzyme level, there is no positive feedback, and hence, no 
bistability. This is manifested in Fig. [TUh by the formation of a cusp at which 
the two thresholds merge. 
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The qualitative variation of the thresholds with the glucose concentration can 
be inferred from Fig. [TOa . Indeed, as the glucose concentration is increased 
from to 1000 /iM, 5 m decreases from S m0 to S moo , shown by the dotted lines 
in the figure. Evidently, as G increases, both thresholds increase, but the ratio 
of the on to off threshold decreases (in contrast to the increase observed in the 
absence of diffusive flux) . 

The quantitative variation of the thresholds with G is obtained by observing 
that (USD yields 



<t>{G) 



fimofy 



1 + 



«1 



G(x) = 4> 



-1 



(36) 



The locus of the thresholds on the Gcr-plane is therefore given by the para- 
metric representation, {G(x), cr(x)}, where G(x) and cr(x) are given by fl36l) 
and (jMj) . respectively. Figure [TUb shows that this representation yields two 
increasing curves that approach each other as G increases. 



3.2.3 Carrier efflux, but no diffusive influx 
In this case, (l2~Ti) becomes 



= S m (G)f(x)a/(K 1 + a) 
X l + 8 m (G)f( X )/K 2 ' 1 ] 

which can be solved for a to obtain 

X[l/QW) + l/«2] ,„ R , 

a x, G) = «i- /aii/ r 38 

The variation of the steady state e and x with a is given by the parametric 
curves, {<r(x, G), e(x, G)} and {a(x, G),x)}, respectively, where e(%, G) and 
a(x, G) are given by (fT7l) and (|38l) . respectively. 

Figure [TTh shows the variation of the enzyme activity with o at G = in 
the presence (full curve) and absence (dashed curve) of the carrier influx. Ev- 
idently, the carrier efflux has no effect on the on threshold, but it significantly 
increases the off threshold. This is because k 2 is so large that the carrier ef- 
flux is negligible in non-induced cells, which contain relatively small enzyme 
and inducer levels (Fig. [Tib). However, the carrier efflux significantly reduces 
the the intracellular TMG level of the induced cells. The net effect of this 
reduction is to increase the off threshold of the induced cells. 

It is instructive once again to study the variation of the thresholds with S m . The 
thresholds satisfy the equilibrium condition ( J37l) and the bifurcation condition 
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(a) (b) 

Figure 11. Comparison of the steady states at G = in the presence (full lines) and 
absence (dashed lines) of the carrier efflux, (a) Variation of the steady state enzyme 
activity with the extracellular TMG level, (b) Variation of the intracellular TMG 
concentration with the extracellular TMG level. 

(122]). which can be solved for 5 m and a to obtain 

Ux) = ^ (4 s " 1 ), (39) 




<40) 

These relations define the parametric representation, {<5 m (x), °"(x)}> °f the 
bifurcation curve on the <5 m <r-plane. 

Figure [T2a compares the thresholds in the presence of carrier efflux (full curves) 
with those obtained in its absence (dashed curves). Evidently, when 5 m < 100, 
the carrier efflux has no effect on both thresholds. Under this condition, the 
enzyme activity and the intracellular TMG level are so small that carrier efflux 
is negligible even in the induced cells. As 5 m increases, carrier efflux from 
the induced cells becomes more significant, and the off threshold becomes 
progressively larger than that predicted in the absence of the carrier flux. 
When 8 m > 10 4 , the enzyme activity is so large that even the non-induced 
cells are subject to carrier efflux. At sufficiently large S m , the off threshold 
merges with the on threshold at a cusp, beyond which there is no bistability. 
The bistability disappears because under this condition, the intracellular TMG 
level is 

a 

X W K 2 ; , 

which is independent of the enzyme level. The destabilizing effect of positive 
feedback therefore vanishes, and the prospect of bistability is eliminated. 

Once again, the qualitative variation of the thresholds with G can be inferred 
from Fig. [T2h . As G increases from to 1000 /iM, 5 m decreases from 5 m0 to 
S moo , shown by the dotted lines in the figure. Evidently, as G increases, so 
do the on and off thresholds, as well as the ratio of the on to off threshold. 
Furthermore, since diffusive influx is absent, there is irreversible bistability: 
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Figure 12. Variation of the on and off thresholds in the presence of the carrier efflux, 
(a) Variation of the on and off thresholds with 5 m . The lower and upper full curves, 
defined by ([39}) and (|30"]) . represent the off and on thresholds, respectively, in the 
presence of the carrier efflux. The full circle shows the cusp at which the thresholds 
merge. The lower and upper dashed curves, defined by (|3ip . represent the off and on 
thresholds, respectively, in the absence of the carrier efflux. The dotted vertical lines 
show the values of 5 m at G = and G = 1000 /iM. (b) Variation of the off threshold 
(lower curve) and on threshold (upper curve) with the glucose concentration, G. 

The off threshold exists for all G, whereas the on threshold disappears at some 
G < 1000 //M. 

The quantitative variation of the thresholds with G follows from ( |39l) . which 
yields 



0(G) 



«2 ( Xfx 



G(x) = F 1 



«2 



Xfx 

fimof I / 



- 1 



(41) 



The variation of the thresholds with G is given by the parametric representa- 
tion, {G(x),cr(x)}, where G(x) and <r(x) ar e given by (|4T|) and (jlO]) . respec- 
tively. This representation defines two increasing curves departing from each 
other (Fig. [T2b). 

We conclude by observing that the geometry of the thresholds (full curves in 
Figure fT2h) does not depend on the precise values of the parameters used in 
the simulations (see Appendix [TJ for details). 



3.2.4 Diffusive influx and carrier efflux 

In this case, the steady state intracellular TMG level is given by fl2~Tl) . which 
can be solved for a to obtain 



<X,G) = - 



- («1 + S m f - p) + V («U + 5 m f - p) + 4Kip 



(42) 



where p = x (1 + $mf ' I ' K 2)- The variation of the steady state e and x with a is 
given by the parametric curves, {cr(x, G), e(x, G)} and {a(x,G),x)}, respec- 
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(a) (b) 

Figure 13. Comparison of the steady states of the general model (full curves) with 
the steady states in the presence of diffusive flux (short-dashed curves) and carrier 
efflux (long-dashed curves), (a) Variation of the enzyme activity with the extra- 
cellular TMG level, (b) Variation of the intracellular TMG concentration with the 
extracellular TMG level. 



tively, where e(x,G) and a(x, G) are given by ffTT|) and (]42]), respectively. 



Fig. [13] compares the steady state profiles obtained at G = in the presence 
of the diffusive influx and carrier efflux (full curves) with those obtained in the 
presence of the diffusive influx only (short-dashed curves) and the carrier efflux 
only (long-dashed curves). It is evident that the upper and lower branches 
of the steady state curve, representing the induced and non-induced cells, 
respectively, are a composite of the steady states obtained in the presence of 
the diffusive influx and carrier efflux. Indeed, the lower branch of the full curve 
coincides with the lower branch of short-dashed curve, and the upper branch 
of the full curve coincides with the upper branch of the long-dashed curve. It 
follows that the steady states and threshold of the non- induced (resp., induced) 
cells is well approximated by the model accounting for only the diffusive flux 
(resp., carrier efflux). The reason for this is clear from the above discussion 
of the special cases. In the absence of glucose, the diffusive influx influences 
the TMG accumulation of non-induced cells only — it is vanishingly small 
in induced cells. Likewise, carrier efflux influences the TMG accumulation of 
induced cells only — it plays no role in non-induced cells. 



The foregoing arguments were based on the steady state profiles obtained 
at G = 0, i.e., 5 m = 5 mQ . It turns out that they are true at almost all 5 m . 
This becomes clear if we compare the thresholds for the general case with the 
thresholds obtained in the presence of diffusive influx or carrier efflux. The 
thresholds for the general case satisfy (121]) and (1221) . which can be rewritten 
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Figure 14. Variation of the on and off thresholds in the presence of the diffusive influx 
and the carrier efflux, (a) Variation of the on and off thresholds with 5 m . The lower 
and upper full curves, defined by (fl3j) and (f4"l"|) . represent the off and on thresholds, 
respectively, in the presence of diffusive influx and carrier efflux. The full circles 
show the cusps at which the thresholds merge. The short-dashed curves, defined by 
(|34h and (|35h . represent the thresholds in the presence of the diffusive influx only. 
The long-dashed curves, defined by (j39H and (|40f) . represent the thresholds in the 
presence of the carrier efflux only. The dotted vertical lines show the values of 5 m 
at G = and G = 1000 /iM. (b) Variation of the on and off thresholds with the 
glucose concentration, G. 
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The first equation can be solved to obtain a(x), which can then be substituted 
in f|44|) to obtain S m (x)- These functions determine the thresholds via the 
parametric representation, {<5 m (x), °"(x)}- 

Figure [l~4"a compares the thresholds thus obtained (full curves) with the thresh- 
olds obtained in the presence of diffusive influx (short-dashed curves) and car- 
rier efflux (long-dashed curves). The thresholds exist only on a finite interval 
- they terminate in a cusp at both small and large 5 m . Over almost its en- 
tire range of existence, the upper branch of the full curve overlaps with the 
upper branch of the short-dashed curve, which represents the on threshold in 
the presence of the diffusive flux. Likewise, the lower branch of the full curve 
coincides with the lower branch of the long-dashed curve, which represents 
the off threshold in the presence of the carrier flux. Thus, the thresholds of 
the general case are a composite of the two special cases. The on (resp., off) 
threshold is well approximated by (|34|) and (|35l) [resp., (|39|) and (HQ]) ], which 
represent the thresholds in the presence of only the diffusive flux (resp., carrier 
efflux) . 



27 



0.2 0.4 0.6 0.8 
Scaled enzyme activity 



Figure 15. Comparison of the observed and predicted effects of inducer exclu- 
sion in S. typhimurium SB 3939 transfected with a plasmid encoding the lac 
operon ([Mitchell et al.l . ll98l Fig. 5A). The open triangles show the steady state in- 
tracellular TMG level in cells induced to various levels, and exposed to a carbon- free 
medium containing 1 mM TMG. The full triangles show the intracellular TMG level 
attained if the medium contains 1 mM TMG and 10 mM a-methylglucoside. The 
full curve shows the fit to the data obtained with Eq. (I45|) and the parameter val- 
ues, 4>s\ G=0 = 1, k x = 0.14 min" 1 , s = 1 mM, K x = 0.8 mM, and K 2 = 41 mM. 
The dashed curve shows the fit to the data with the same equation and parameter 
values, the only difference being that <t ) s\G=iomM = 0-5. 

It is evident from Fig. [T4"k that as G increases, both thresholds increase, but the 
ratio of on to off thresholds can increase, decrease, or pass through a maximum. 
For the parameter values shown in Table CD, 5 m o and 6 moo are such that the 
ratio passes through a maximum. The precise variation of the thresholds with 
G is given by the parametric curve, where 



1 



G( x ) = r 1 



a relation that follows immediately from (j44]) . This parametric curve consists of 
two increasing curves that diverge at low glucose concentrations, and converge 
at high glucose concentrations (Fig. H4b). 



3.2.5 Comparison of the model simulations with the data for TMG 

Comparison of the model simulations (dashed lines in Fig. [2h>) with the data 
shows that the simulated off thresholds are higher, and the simulated on 
thresholds are lower, than the corresponding thresholds observed experimen- 
tally. Part of this discrepancy probably stems from the fact that the parameter 
values, which correspond to various strains of E. coli, are somewhat different 
from those shown in Table [Q However, the analysis of the data suggests that 
there are two additional sources of the discrepancy. 



Dependence of S on the the enzyme activity The model assumes that 
the effect of inducer exclusion, characterized by the function <p s , is completely 
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determined by the extracellular glucose level. However, the data in Fig. 
implies that <p s also depends on the activity of the permease. The open tri- 
angles in the figure show the steady state intracellular TMG levels in cells 
induced to various levels, and then exposed to a carbon-free medium contain- 
ing 1 niM TMG. The full triangles show the steady state intracellular levels 
attained when the carbon-free medium contains 1 mM TMG and 10 mM a- 
methylglucoside (aMG), a non-metabolizable analog of glucose that mimics 
its inducer exclusion effect. It is clear from the data that inducer exclusion 
is significant in cells containing low enzyme levels, but it disappears in fully 
induced cells. 

If cf) s is completely determined by the extracellular glucose concentration, there 
is no combination of model parameters that can capture the data in Fig. [15j To 
see this, observe that the data was obtained in a carbon-free medium, which 
prevents efficient enzyme synthesis. The appropriate equation for modeling 
the data is therefore given by the following variant of Eq. ([5]) 

Since the extracellular TMG level, s, is constant (1 mM), the above equation 
implies that at any given x, the enzyme activity in the absence of aMG (or 
glucose) is 4> S {G) times the activity in the presence of aMG. This contradicts 
the data. At x ~ 30 mM, for instance, the enzyme activity in the absence of 
aMG is roughly half the activity in the presence of 10 mM aMG. However, at 
x ~ 40 mM, the activities are the same in the absence and presence of aMG. 

The molecular explanation for the data in Fig. [15] sug gests that 4> K depend s 



not only on G, but also on the prevailing enzyme level (jMitchell et all Il982l ). 
Indeed, inducer exclusion occurs because the uptake of aMG or glucose results 
in the formation of dephosphorylated enzyme IIA glc (a transport enzyme for 
glucose), which inhibits lac permease by binding to it. In fully induced cells, 
the permease level is so high that most of it is essentially free of IIA glc . As the 
permease level decreases, so does the inhibition due to LacY-IIA glc bindingO] 
Thus, the magnitude of the inducer exclusion effect depends on the prevailing 
permease level. 

The model, which ignores the dependence of S on the permease level, pro- 
vides a good fit to the data obtained in the absence of a-MG (full line in 
Fig. [T5|) . However, it fails to capture the data obtained in the presence of a- 
MG (dashed line in Fig. [T5|) . It overestimates the inducer exclusion effect in 
induced cells, and underestimates this effect in cells containing lower enzyme 
levels. Consequently, the observed off threshold is expected to be lower, and 



3 No data was obtained at scaled enzyme activities below ~ 0.25 because the re- 
combinant strain used in these experiments was quasi-constitutive. 
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Figure 16. Variation of the specific growth rate due to induction of lac. (a) During 
growth of E. coli B in the presence o f succinate and TMG, full induction of the cells 
reduces the specific growth by ~30% (jNovick and Weinerl . ll957l . Table 2). (b) During 
growth of E. coli K12 MG1655 in the presence of glycerol and IPTG, full induction 
reduces the specific growth rate by only ~5% (iDekel and Alonl . 120051 . Fig. 2a). 



the observed on threshold is expected to be higher, than the corresponding 
thresholds predicted by the model. 



Dependence of r g0 on the enzyme activity The model assumes that 
in the absence of glucose, the specific growth rate is a fixed constant, r 9i0 - 
However, the data shows that the specific growth rate varies significantly with 
the enzyme level of the cells (Fig. [T6h). When the cells are exposed to 0.5 mM 
TMG, the specific growth rate is ~30% lower than that of non-induced cells, 
and the graph suggests that the specific growth rate declines further at higher 
extracellular TMG levels. The reduced specific growth rate of the induced cells 
serves to decrease the off threshold because the lower the specific growth rate, 
the higher the enzyme and intracellular TMG levels. 



In ea rlier work, we studied the dynamics of growth on lactose (INarang and Pilyuginl . 



20071 ). There, we showed that bistability is suppressed during growth on lac- 
tose because the specific growth rate increases with the enzyme level, thus 
enhancing the stabilizing effect of dilution. Since the specific growth rate de- 
creases with the enzyme level during growth on TMG, it is plausible to expect 
that in this case, the stabilizing effect of dilution is depressed, and the bistable 
regime is enlarged. 



3.2.6 Comparison of the model simulations with the data for IPTG 

The above arguments suggest that the discrepancy between the data and the 
simulations would be reduced if the experiment was performed in the absence 
of glucose with a gratuitous inducer that does not result in a significant reduc- 
tion of the specific growth rate. This conclusion seems to be consistent with 
the data. 
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Figure 17. Comparison of the data and the simulation for bistability with IPTG. 
(a) Variation of the observed steady sta te /3-galactosidase a ctivity of E. coli K12 



with extracellular IPTG levels at G = ([Laurent et all 120051 . Fig. 7) . (b) Variation 



of the simulated enzyme activity with extracellular IPTG levels at G = 0. The 
dissociation constant, K% , was assumed to be 7 /j,M. All other parameters were 
assumed to be the same as those used for simulating the experiments with TMG. 

It turns out that when the cells are grown in the presence of IPTG, the specific 
growth rate of the induced cells is only ~5% smaller than the specific growth 
rate of the non- induced cells (Fig. [T6b). Recently, Laurent et al. measured the 
variation of the enzyme activity with the extracellular IPTG levels in the ab- 
sence of glucose (Fig.fTTk). They found that the mean activity does not display 
the discrete all-or-none response observed by Ozbudak et al. (presumably due 
to the physiological heterogeneity of the cells and the inherent stochasticity 
of the induction process). However, there is clear evidence of bistability in 
the 5-fold range, 0.3-1.5 /iM, which is significantly smaller than the 10-fold 
range shown in Fig. [2|l. The IPTG concentrations supporting bistability are 
significantly smaller than those observed by Ozbudak et al — the off and on 
thresholds are 0.3 and 1.5 fiM (compared to 3 and 30 /iM in Fig. [2^l) - This 
is because the affinity of IPTG for the repressor is an order of magnitude 
higher than that for TMG. Simulation of the model with K~ x = 7 yields 
results that are in reasonable agreement with the data obtained by Laurent et 
al (Fig.Eb). 

It remains to be seen if an extension of the model, accounting for the depen- 
dence of (f) g and r 9i o on the prevailing enzyme level, can resolve the discrepancy 
with respect to the data for TMG. 



4 Conclusions 



Analysis of the experimental data shows that the diffusive influx and carrier 
efflux have a profound effect on inducer accumulation in non-induced and 
induced cells, respectively. Since bistability entails the coexistence of steady 
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states corresponding to both non-induced and induced cells, neither one of 
these fluxes can be neglected. We analyzed a model of lac bistability taking 
due account of the diffusive influx and carrier efflux. We find that: 

(1) The diffusive influx decreases the on threshold. 

(2) The carrier influx increases the off threshold. 

(3) Over a wide range of permease activities, the diffusive influx has no ef- 
fect on the off threshold, and the carrier efflux has no effect on the on 
threshold. Since each of these fluxes influences only one of the thresholds, 
the combined effect of these fluxes can be captured by simple analytical 
expressions corresponding to two limiting cases. 

(4) Simulations of the model show good agreement with the data for IPTG. 

(5) There is a significant discrepancy with respect to the data for TMG. 
Analysis of the data suggests that this discrepancy is due to the depen- 
dence of the inducer exclusion effect and the specific growth rate on the 
prevailing lactose enzyme levels. 
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A Necessary condition for a bifurcation 



The Jacobian of Eqs. fflll) - (fT4|) at any steady state is 
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Evidently, two of the eigenvalues are — <j) g . Hence, the steady state is stable if 
and only if the remaining two eigenvalues corresponding to the submatrix 



have negative real parts. Since the trace is always negative, stability results 
whenever the determinant is positive. Hence, the necessary condition for a 
bifurcation is that the determinant be zero, that is, 
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Figure B.l. There are exactly two lines passing through the origin (dashed lines) 
that are tangential to the graph of / (full curve). The slopes of these lines are 
f(xi)/xi = fxiXi) and f{x2)/X2 = f x (X2)- 

Since (j) 2 g <C 4> 9 /t x , the bifurcation condition can be approximated by the 
relation 
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where = s e /0 3 and S m = 05 m o- This is the same as Eq. ( !22j) . 

The asymptotic dynamics of Eqs. (ITT|) - (fHj) reduces to the stability analysis of 
its steady states. Indeed, Eqs. (fTTi) - (fl2i) are uncoupled from equations (fTBl - 
( JT4l) . A direct calculation shows that the subsystem (TIT]) - (fl2l) has a negative 
divergence (i.e. the trace of the Jacobian is negative) in the entire state space. 
An application of the Bendixson-Dulac criterion eliminates the possibility of 
periodic orbits and saddle connections, hence any solution (x, e) of (TTTT) — (TX2T) 
converges to some steady state. It follows immediately that the solution (7, p) 
of ffT3]) - ffl"4|) also converges to the corresponding steady state. 



B The positive roots of Eq. (J30D 



We wish to show that q(x) = X ~ f / fx nas exactly 2 positive roots, \i < X2, 
and is positive if and only if xi < X < Xz- This assertion has the a simple 
geometric interpretation. Since f(x) is a sigmoidal function for sufficiently 
large a and a, there are only two lines passing through the origin that touch 
the graph of /, and they do so at the two points, xi an d Xi (Fig. IB.lj) . 
Furthermore, for any Xi < X < X2, the slope of the line passing through (0, 0) 
and (x, fix)) is smaller than f x . 

The formal proof is obtained by observing that q(0) = — /(0)// x (0) < 0, and 
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Figure C.l. Geometry of the S m (x) and cr(x) in the case of diffusive influx, 
at large x, 
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It follows that has at least 2 positive roots. In fact, it has exactly 2 
positive roots because it has only 1 extremum on [0, oo). Indeed, 



dq = 
dx 



fxxf fy _ fxxf 



J Y 



J Y 



which is zero if and only if f xx = 0, i.e., x — X*i where (x*^f(x*)) is the 
unique inflection point of the sigmoidal function, /. 



C Bifurcation diagram for diffusive influx 



We wish to show that regardless of the parameter values, the functions, cr(x) 
and S m (x), defined by (|34j) and (|35l) . respectively, have the geometry shown 
in Fig. IC.ll It follows that as X increases from X i-> ° increases and 5 m de- 
creases until a cusp forms at the point X = X* where da/dx = dS m /dx = 0. 
Beyond x = X*i a decreases and S m increases until they approach and oo, 
respectively. Thus, the variation of the thresholds on the 5 m a-plane has the 
geometry shown in Fig. [TUh . regardless of the parameter values. 

Now, the geometry of u{x) is identical to that of q(x), which was analyzed in 
Appendix [B] Hence, it suffices to show that S m (x) has the geometry shown 
in Fig. IC.ll To this end, observe that (|35|) implies \im x ^ x + x - 5 m (x) — °°- 
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Figure D.l. Geometry of the 6 m (x) and <r(x) in the case of carrier efflux. 
Moreover, since 
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dS m /dx and da/dx have opposite signs at all Xi < X < X2, except when 
X = x*i at which point dS m /dx = da/dx = 0. 



D Bifurcation diagram for carrier influx 



We show below that regardless of the parameter values, S m (x) and a(x), as 
defined by (139|) and (flOl) have the geometry shown in Fig. ID.ll Hence, the 
variation of the thresholds on the 5 m cr-plane always has the geometry shown 
in Fig. r 



To see this, observe that since q(x) = X ~ fx/ f nas exactly two positive 
zeros, xi < Xi (Appendix [B]), so does S m (x)- Furthermore, S m (x) has a unique 
maximum on xi < X < Xi- Indeed, if we let z = x + 1, ^m(x) can be rewritten 

as 
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which implies that 
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Now, the polynomial b(z) changes sign exactly once in the interval (1, +00) 
because 6(1) = 4a + 8a > 0, b(z) < for large z, and b zz (z) = 6a (1 — z) < 
for all z > 1 (i.e. b(z) is a concave function). Thus, the function S m (x) has a 
unique maximum. 

The function, c(x)> becomes infinite at points, < xi < Xu, lying within 
(Xi,X2)- Indeed, fl40]) yields 

^1 = fornix) _ 1 

so that c — > 00 if and only if 

X 2 f x .. fxfx A Xf x f, X 



l^-^Vf'j)^ f V K 2 

Now, S m (x) becomes zero if and only if xfx/ f = 1- ^ follows that k 2 is 
sufficiently large, a — > 00 at two points, X2 < X« lying within (xi, X%)- 

Finally, cr(x) has exactly 1 extremum in (xu Xu) which is attained at the very 
same value of x that maximizes 5 m . Indeed, 



Ki da ( a x \ dS. 

0. 



in 



m (K 1 + a) 2 dx XKx + a k 2 J dx' 
where the term in parentheses is positive whenever S m (x) > because 

g(x) = x 2 f g(x) _ X_ = xf 

Ki + a(x) K 2 (xf x ~ f) K i + a {x) k 2 K 2 (xf x ~ f) ' 

It follows that d5 m /dx and da/dx have opposite signs in the interval Xi < 
X < X21 and the functions 5 m and a attain their extrema at the same values 
of x- Since S m attains a unique maximum in the interval xi < X < X2, so does 
a. 
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